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We demonstrate for the first time natural phase matching for optical frequency doubling in a 
high-Q whispering gallery mode resonator made of Lithium Niobate. A conversion efficiency of 
9% is achieved at 30 fj,W in-coupled continuous wave pump power. The observed saturation pump 
power of 3.2 mW is almost two orders of magnitude lower than the state-of-the-art. This suggests 
an application of our frequency doubler as a source of non-classical light requiring only a low-power 
pump, which easily can be quantum noise limited. Our theoretical analysis of the three-wave mixing 
in a whispering gallery mode resonator provides the relative conversion efficiencies for frequency 
doubling in various modes. 

PACS numbers: 42.65.Ky, 42.60.Da 
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Various types of oscillations can be enhanced by ap- 
propriate resonators. Examples range from resonating 
chambers for sound waves to cavities for electromagnetic 
waves. Nonlinear cross-coupling of different oscillations 
can be enhanced as well, e.g. between optical waves, be- 
tween optical and radio- frequency electromagnetic waves, 
and, a subject of considerable recent attention, between 
mechanical and optical oscillations [Ij. 

A resonator of particular interest in the context of 
quantum and nonlinear optics is the whispering gallery 
mode resonator (WGMR). In WGMRs, light is guided 
via continuous total internal reflection. The small mode 
volume combined with the absence of external reflecting 
elements provides for a very high quality factor Q and 
good spatial localization of the optical fields. In addi- 
tion to these important capabilities, the WGMR enables 
us to continuously control the coupling strength between 
the intra-cavity modes and the modes outside the cavity. 
This capability may be compared to a cavity mirror with 
continously variable reflectivity. 

Due to these unique capabilities, various nonlinear pro- 
cesses have been observed in WGMRs, including four 
wave mixing jll-Q , Raman [sl, Q , parametric j? [and Bril- 
louin Q scattering, microwave upconversion 3], second 
[lOj and third [llj harmonic generation. Furthermore, a 
plethora of dynamical effects has been predicted in the 
cavities with the second order nonlinear response, exam- 
ples being bistability, chaotic behavior and self pulsing 
[13] , along with non classical effects such as squeezing and 
entanglement [13| . It is intriguing to investigate these ef- 
fects in a highly efficient WGMR. 

In this paper we report on the second harmonic (SH) 
generation in a doubly resonant WGMR. We employ non- 
critical type-I phase matching, and achieve the SH con- 
version efficiency of 9% at 30 fiW pump power. This is 
nearly two orders of magnitude higher than the previous 
state-of-the-art result for the same pump power, obtained 
with a less efficient quasi phase matching fl^ . At higher 



pump power, the conversion efficiency saturates due to a 
dynamic "self-limiting" effect. 

The phase matching conditions for SH and degener- 
ate parametric down conversion (PDG) reflect the con- 
servation laws associated with the resonator symmetries. 
They arise from the interaction Hamiltonian 



Hi, 



J x^^^E,{ElfdV + h.c. 



(1) 



expressed via the effective quadratic nonlinearity x^^^ 
and the pump and SH electric fields Ep and Eg. In 
Cartesian coordinates, the integral in ([Ij gives rise to 
three momentum conservation equations: 2kp — kg. In 
a spherical WGMR, a mode is defined by a set of po- 
lar, azimuthal and radial numbers {L, m, g}. Assuming 
that the WGMR radius R is much larger than the optical 
wavelength A, such that there is no coupling disturbance, 
and that the WGMs are near-equatorial, we can reduce 
the electric fields in both polarizations to a scalar eigen- 
function [lj| with a scaling factor Eq : 



ELmqir, 9, 4,) = EqYlUO, 4)jL{kLqr). (2) 
The angular part of ([2]) is a spherical harmonic 
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where P™(cos6') are the associated Legendre polynomi- 
als. The radial wave number k^q in the argument of 
the spherical Bessel function jL{kLqr) relates the WGM 
eigenfrequency a; to L and q via the dispersion equation 



L r 
R 
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where n is the material index of refraction and aq is the 
q-th zero of the Airy function [l3, [iBl • 



FIG. 1: (Color online) Normalized Clebsch- 

Gordan coefficients (m + Up,m + Vp;m,m\2m + 
Vs, 2m)/(m, m; m, m\2m, 2m). 

Applying ((T|) to the angular part of ^ yields the 
Clebsch-Gordan coefficients whose subscripts s and p re- 
fer to the pump and the second harmonic, respectively: 

{Lp, Lp;mp,mp \ Ls,ms) ^ (5) 

= / Yl.r. id, cj,)Yl,,^ {9, 4>) sin OdOdcj,. 

Selection rules for Clebsch-Gordan coefficients are well- 
known. In our case: 

nig — 2mp, Lg < 2Lp. (6) 

The first condition in ^ is conservation of the azimuthal 
component of the orbital momentum, while the second 
one arises from the triangular condition for a vector sum 
of the orbital momenta of two pump photons. The phase 
matching conditions for the SH generated on a surface of 
a spherical WGMR have been already construed as the 
addition rules for the orbital momenta . This also ap- 
plies to our WGMR, and the Clebsch-Gordan factors ([5]) 
in the conversion efficiency emphasize the analogy with 
atomic systems. Notice that a spectrum of a spherical 
WGMR (HI) is degenerate with respect to to, which de- 
couples the first selection rule of (O from the frequency- 
doubling resonance condition 

ixjg = 2ujp (7) 

and makes the phase matching easier to satisfy. 

The relative SH conversion efficiencies for different 
WGMs are found from Eq. ((5)) using Gaunt's formula [l^ 
and treating large factorials in Stirling's approximation. 
Remarkably, for L and to sufficiently large for dispersion 
equation Q to remain a good approximation (which is 
the case in our experiment), the coefficients ([5|) normal- 
ized to (to, to; to, to|2to, 2to) are practically independent 
of TO and therefore of R and A. This normalization co- 
efficient describes nonlinear coupling of two equatorial 



FIG. 2: (Color online) Normalized absolute-value coupling 
coefficients for frequency doubling {ajp,gp} — >■ {2ijjp,qs}. 

modes {L = to). The equatorial modes have a single anti- 
node at = 7r/2 and therefore yield the best coupling of 
the WGMR with the free-space modes, as well as the 
best overlap with other equatorial modes. To describe 
non-equatorial modes, we introduce Vp = Lp — rup > 
and Vg = Lg — nig > 0. The normalized Clebsch-Gordan 
coefficients are shown in Fig. [TJ From this Figure we 
see that the second selection rule of ^ breaks the sym- 
metry between the SH and degenerate PDC processes. 
While the former can couple an equatorial pump only to 
another equatorial mode, the latter can couple the equa- 
torial pump to a variety of modes. 

The radial part of integral ([1]) with substituted ([2]) 
and (HI) is evaluated approximating the large-order Bessel 
functions by Airy functions. The results are normalized 
to the {ujp,qp = 1} — 5- {2ujp,qs — 1} coupling coeffi- 
cient, and shown in Fig. [2j These coefficients are practi- 
cally independent of Lp and Lg and therefore of R and 
A. The relative efficiencies of different conversion chan- 
nels {cjp, Lp, TOp, qp} — > {2u!p, Lg, 2TOp, qg} are found as a 
product of squares of coefficients from Figs. [T]and[2] 

The multitude of conversion channels allows for a con- 
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FIG. 3: (Color online) Experimental setup. 
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FIG. 4: (Color online) Transmitted pump power and the SH 
signal. The coupling is optimized for high SH output, which 
corresponds to over-coupled fundamental modes. On the in- 
set: a typical WGM spectrum of a critically coupled resonator 
at the pump wavelength. 

venient experimental investigation. Our experimental 
setup is shown in Fig. [31 The WGMR is made from a 
congruent 5% MgO-doped Lithium Niobate z-cut wafer. 
Its radius is i? = 1.9 mm and height is approximately 0.5 
mm. The rim of the disk is polished into a tapered shape 
that near the equator can be approximated as a spheroid 
with radii R and p. Their ratio R/ p w 7.5 is nearly opti- 
mal for coupling of the horizontally polarized pump. 

The pump is a free-space Gaussian beam from a con- 
tinuous wave Nd:YAG laser (InnoLight Mephisto) with 
a wavelength of 1064 nm. It is focused by a microscope 
objective onto the back side of a diamond prism at an 
angle that allows for total internal reflection. The res- 
onator is placed near the reflected beam's footprint, so 
that the in- and out-coupling occurs via frustrated total 
internal reflection. The evanescent field coupling is varied 
by changing the gap between the prism and the resonator 
rim. The disk and the prism are mounted on a hot plate 
whose temperature is stabilized. The horizontally polar- 
ized pump couples to the ordinary polarized modes of the 
resonator, and the vertically polarized SH is generated 
into an extraordinary polarized mode. Light coupled out 
of the disk is re-collimated by another objective, and the 
fundamental and the SH components are separated by a 
dichroic mirror and focused onto photodetectors. Finally, 
the signals are monitored by an oscilloscope. 

Scanning the laser frequency we observe the pump 
transmission, which reveals the WGM spectrum at the 
fundamental wavelength (see Fig. Fitting the res- 
onances with a Lorentzian function we find a minimal 
linewidth 7 = 8.3 MHz, (Q = 3.4 x lO'^), and a maxi- 
mum contrast of 90%, for a critically coupled resonator. 

An electrode attached to the disk allows us to apply a 
bias voltage to the resonator. The bias shifts the ordi- 
nary (pump) and extraordinary (signal) mode frequencies 
Up and ujg differently, proportionally to the electro optic 
coefficients r^i and r^s w 3.2r3i respectively. This dif- 
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FIG. 5: (Color online) SH conversion efficiency as a function 
of the in-coupled power: experimental data and theoretical fit 
[lol ] corresponding to the saturation pump power Wo — 3.2 
mW. On the inset: transmitted pump power and the sec- 
ond harmonic signal above the "self-limiting" threshold which 
arises at a much lower power than Wq. 

ferential frequency-tuning of WGMs allows us to achieve 
condition ([7]) required for the frequency doubling. 

Changing the resonator temperature has the same ef- 
fect as changing the bias voltage, since the temperature 
dependences of the indices of refraction no{T) and iie{T) 
are different. While the electro-optical tuning is faster, 
the temperature tuning has a larger dynamic range. 

To find the phase matching temperature for the 
WGMR, we made a bulk crystal from the same wafer 
and determined its phase matching temperature to be 
121.7°C. This result was consistent with the calculations 
based on the Sellmeier equation jl7j and with an ear- 
lier experiment Ts*]. Using the Sellmeier equation to- 
gether with the spherical WGMR dispersion 1^, we cal- 
culated the optimal phase matching temperature to be 
94.1°C. We assumed the following parameters: ujp = 
27rc/1064.3nm, nip « 25200, and Qs = Qp = 1, as this 
is the most efficient conversion channel within the acces- 
sible temperature and wavelength ranges. 

We find the phase matched modes by frequency tun- 
ing the extraordinary WGMs at 532 nm relative to the 
ordinary WGMs at 1064 nm with the bias voltage. We 
compensate for the frequency shift of the fundamental 
WGM by the laser central frequency. At some bias volt- 
ages the requirement ([7]) is occasionally satisfied for a 
pair of modes. Then if the phase- matching conditions 
([6]) are also satisfied, SH is generated, see Fig. |4l For 
different pump modes this might happen once or several 
times. However for the majority of the modes it never 
happens within the explored parameter space. 

We observed different conversion efficiency for differ- 
ent WGMs. Fig. [5] shows the best measured conversion 
efficiency as a function of the in-coupled pump power. 
It is worth noting that the detected SH power is propor- 
tional to the power circulating in the resonator and to the 
resonator coupling. The circulating power, on the other 
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hand, is reduced when the couphng increases (due to de- 
creasing Q- factors). As a result, the detected SH power 
as a function of the couphng has a maximum which cor- 
responds to the over-couplcd fundamental WGM. 

The theoretical model for the SH generation [19J in 
an idealized doubly resonant WGMR has been advanced 
in [lo| . A key parameter in this model is the satura- 
tion power Wq which can be considered as an absolute 
measure of the nonlinear conversion efficiency, indepen- 
dent of the pump power. The experimental value of 
Wq reported in is 300 mW. In our experiment we 
measured Wq = 3.2 mW, see Fig. [S] However even 
this low power greatly exceeds the theoretical estimate 
of Wo !=s 6 fiW derived for our resonator based on the 
theoretical model in rcf. fio]. This indicates that the 
spatial overlap between the SH and fundamental WGMs 
in our experiment deviated from the optimal one by a 
factor of ^yWo {theo)/Wo (exp) w 0.043. To explain this 
discrepancy we point out that the diagonal coupling co- 
efficients Lg = 2Lp in Fig. [T] decrease rather slowly with 
increasing L — m. Also, according to Fig. ^ coupling 
between WGMs with different radial numbers q remains 
significant even when the difference is large. Therefore we 
observe multiple instances of phase matching with non- 
optimal conversion efficiencies and consequently higher 
saturation power Wq, than for optimal conditions. It is 
difficult to uniquely establish which pair of modes yielded 
the overlap factor inferred from experiment. 

Similar non-optimal conversion efficiencies were mea- 
sured in the range of 94°C to 122°C. Evidently, the cal- 
culations based on the spherical WGMR dispersion 
have not been sufficiently accurate for our non-spherical 
resonator. Gonsidering a very narrow temperature width 
of a WGMR phase matching, of the order of 0.001°C, em- 
pirically searching for the optimal phase matching within 
a 30-degree temperature range was cumbersome. A more 
accurate initial temperature estimate is preferable. 

The conversion efficiencies measured both in [13] and 
here (see Fig. [5]) reach the maxima and saturate at much 
lower in-coupled pump powers than Wq. We attribute 



this to a nonlinear "self-limiting" effect arising when the 
circulating pump power exceeds a certain threshold, see 
the inset in Fig. [5l This effect is not related to the pho- 
torefractive damage (which can also be observed in our 
resonator, see ^]). While the photorefractive damage 
causes quasi-permanent changes to the WGM spectra, 
the "self-limiting" effect shows purely dynamic behavior 
as a function of the pump power and coupling. This ef- 
fect may be related to the onset of self-pulsing or similar 
dynamics [l2| and will be discussed elsewhere. 

In summary, we present the first natural non-critical 
type-I phase matching for 1064 to 532 nm optical fre- 
quency doubling in a high-Q WGMR. This results in a 
hundred fold lower saturation pump power. Our theoret- 
ical study of phase matching in the spherical geometry 
of WGMRs has predicted the existence of multiple cou- 
pled WGM pairs with varying conversion efficiency con- 
sistent with our observations. The study also predicted 
the asymmetry between the SH and PDC selection rules, 
which may lead to the competition between these pro- 
cesses and explain the observed "self-limiting" effect. 

The best efficiency of 9% was achieved at 30 /iW in- 
coupled CW pump power. Such a high-efficiency low- 
power frequency conversion process is very interesting 
in the context of nonlinear opto-mechanical coupling [l| , 
nonlinear dynamics [l2| and quantum optics [13]. Exper- 
iments with linear cavities showed reduction of amplitude 
noise below the quantum limit at the fundamental [2l| 
and SH [3] frequencies. The theoretical limit for this 
reduction has been shown [l^] to be determined by the 
pump noise. In the low pump power regime, shot noise 
and sub-shot noise pump light fields are easily accessible, 
see e.g. [2^. We plan to exploit the low-power opera- 
tion regime in nonlinear WGMRs to lift the pump noise 
limitation and to realize a source of non-classical light. 
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